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1. Introduction 

Soliton equations has been largely studied in the past in terms of integrability conditions 
and solutions |U|2l[3lllJ|5l|6l[7J[8l|9l HOI ED]- These solutions have been discovered 
through numerous techniques such as Backlund transformations, inverse scattering 
methods, Lax pairs formulation, G'/G-expansions, Wronskian formulation and bilinear 
Hirota method. The later ones will be of interest in this paper. To this instance, the 
Hirota derivative is defined as 

K(f ■ 9) = (d Xl - d X2 ) n f(x 1 )g(x 2 )\ x=Xl=X2 . (1) 

In this paper, we consider the Boiti-Leon-Manna-Pempinelli (BLMP) equation 

0, (2) 

where u = u(x,y,t) and subscripts represent partial differentiation with respect to the 
given variable. This equation has received small attention [TT] and here we propose to 
give some new classical solutions through the Hirota and Wronskian methods. To do so, 
we note that the BLMP equation is strongly related to the Korteweg-de Vries (KdV) 
equation. Indeed, let us take 

u(x,y,t) =p(z,t) = p(x + q(y),t), (3) 

where q is an arbitrary function of the variable y, then p satisfies 

Pzt +Pzzzz - QPzPzz = 0. 

Defining h = p z , we see that h satisfies the KdV equation, which has the following 
bilinear form [H [9] 

V z (V t + V 3 z )(f-f) = 0, (4) 

using the change of variable h = —2d z log /. Thus p = —2d z log / satisfies the BLMP 
equation for f(z,t) a solution of the above bilinear equation associated to the KdV 
equation. This suggests that the BLMP equation possesses multisoliton and rational 
similarity solutions and many other solutions as we will demonstrate. 

Without doing the reduction ([3l) in u(x, y, t) but using the change of variable given 
above, we can bilinearize, in general, the BLMP equation (T5]) by casting 

u(x, y, t) = -2d x log f(x, y, t) - m(y), (5) 

where m = m(y) is an arbitrary function of y. One then gets the bilinear equation 

(V y (V t + Vl) + 3m'(y)Vl)(f ■ f) = 0. (6) 

After considering the classical analysis of the BLMP equation, it is natural to ask 
what would be its supersymmetric (SUSY) counterpart. SUSY extensions of nonlinear 
partial differential equation (PDE) has, in the past few decades, been largely studied 
in terms of integrability and solutions, see for example [12], [131 CLU HSJ US, [TTJ, UHl EIHJ 
[2Ql [2H [221 [231 12H I2ZI [211 [291 [30]. Here, we consider for the first time a M = 1 SUSY 
extension of the BLMP equation (T5]). 
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To produce a M = 1 SUSY extension of the BLMP equation, we extend, as usual, 
the space (x,y,t) to the superspace (x,y,t;9) [25] where 9 is a Grassmann variable 
which satisfies 9 2 = 0. As usual, we consider a fermionic superfield y,t;9) to get 
non trivial SUSY extensions [121 [131 [151 HI HE] . A A/" = 1 SUSY extension of the BLMP 
equation (|2]) is given as 

Dfa + D 2 y Dl<$> - 3D x D y (D 3 x $D y D x $) = 0, (7) 
where D x and D y denotes the supercovariant derivatives defined as 

D x = d + 9d x , D y = de + 9d y , D 2 x = d x , D 2 y = d y . (8) 

Using the finite Taylor expansion 

y, t; 9) = £(x, y, t) + 9u(x, y, t), (9) 

where £ is a complex- valued fermionic (odd) function and u is a complex-valued bosonic 
(even) function, we get a set of two equations 

Uyt ~\~ Uxxxy 3u xx lly 3u x U X y -\- 3^ x ^ xxx 0, (10) 
^yt ^xxxy 3u X y^x 3ll x ^ X y 0. (H) 

In the fermionic limit where £ is set to zero, we retrieve the classical BLMP equation as 
expected. 

As in the classical case, the SUSY BLMP equation is related to the SUSY KdV 
equation. Indeed, supposing that £(x,y,t) = \J q'{y)x( z i t) and u(x,y,t) = v(z,t) where 
z = x + q(y), then the above set of equations leads to 

v z t + v zzzz -6v z v zz + 3xzXzzz = 0, (12) 
Xt + Xzzz - 3xzVz = 0. (13) 

In the fermionic limit, we get the classical KdV equation for the variable v z . The 
equations (fT2j) and (fTBl becomes the SUSY KdV equation first proposed by Manin and 
Radul [31] after identifying Xz = V 9 - We thus get a series of solutions of the SUSY 
BLMP equation using this procedure. 

In section 2, we give new classical solutions of the BLMP equation using the 
Wronskian and Hirota bilinear formulation. In section 3, we introduce the SUSY Bell 
polynomials which will be the key ingredient of our analysis. In section 4, we give 
solutions of the SUSY BLMP equation, e.g. soliton and periodic solutions. There is 
also an attempt to find the superpartner solutions of the given classical solutions u of 
the BLMP equation. In section 5, we produce bilinear Backlund transformations using 
the SUSY Bell polynomials described in section 3. 



2. New classical solutions 



In the introduction, we have shown a link between the BLMP equation and the KdV 
equation. It means that, for example, we can produce rational similarity solutions for 
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the BLMP equation based on the corresponding ones for KdV. Indeed, for the rational 
similarity solutions [H [21 111 13 [3 [10] , we get 

u n (x, y, t) = <p n (z, t) = -2— log Q n (z(3t)~^), 

where the functions Q n are known as the Yablonskii-Vorob'ev polynomials and are 
defined recursively [7] as 

Qn+lQn-l = ZQI - 4(Q n Qn,zz ~ Q',J, 

with Qo(z) = 1 and Qi(z) = z. For example, we get the solutions 
, N 6z 2 . , 12(30tz 2 + z 5 ) 

U 2\Zit) = -—— -, U 3 {z,t) - 



12t + z 3 ' y ' ' 720t 2 - 60tz 3 - z 6 
This infinite series of rational similarity solutions have not been obtained before. 

We can exhibit soliton solutions. Indeed, for the iV-soliton solution one gets for the 
function / [HE] 

/at = ex P I + X! > ( 14 ) 

/U=0,1 \i=l i<j / 

where ^ = k^z — n\t and e a!i = Ay = . The above solutions has been deduced 

from the reduction of the BLMP equation to the KdV equation via (J3j). 

Using a Wronskian formulation for the KdV equation [21 El [8], we get rational, 
positon, negaton and complexiton solutions of the BLMP equation. To this instance, 
we look for solutions / = f(z,t) of the following form 

f = det(H^,H^,---,H^), (15) 

where H = = (hi, h 2 , ■ ■ • , h N ) T and H® = for i = 1, 2, • • • , N - 1. The 

resulting function u = —2d z log / is called a Wronskian solution of order N. The 
Wronskian solution requires pi El El E] that the functions hi satisfy the following linear 
PDE's 

— hi t zz — \hi, hij = —4hi^ zzz , (16) 

for i = 1,2, • • -N. When Aj is real, we get three cases of solutions: rational (Aj = 0), 
positon (Aj > 0) and negaton (Aj < 0). For the complexiton solution, the eigenvalues Aj 
are supposed to be complex. 

Here we give some examples and refer the interested readers to [31 El El [11] for 
further discussions on the Wronskian solutions. 

A negaton solution of second order 

ii U t 7l ) = 8 7l cosh 2 ( 7l (.-4 7l 2 t)) 

n{ ' ' 7lj 2 7l (12 7l 2 i-z) -sinh(2 7l (^-4 7l 2 t))' 1 } 

a positon of second order 

u ( Z t 7l ) = 8 7l cos 2 ( 7l (, + 4 7l 2 t)) 

p[ ' ' 7lJ ~ 2 7l (12 7 2 £ + z) + sin(2 7l (^ + 4 7l 2 t)) U8J 
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Figure 1. From left to right, a 2-soliton solution ([H)) for N = 2 (k\ = \, k% = 1) 
and t = 20, a negaton of second order (|17p (71 = 1) and a positon of second order (fT5|) 
(71 = 1) for 2 = x + q(y) and t = 2. 



and a complexiton 

u c (z,t,r]) 



4ilm(i]) cosh(x) cosh(%) 



(19) 



v /r]cosh(x) sinh(x) — y/rj cosh(x) si nn (x) 

where x = \/v{ z ~~ 4r/t), x = V 7 ^- 2 — ^t) and rj = a + i/3. There is numerously more 
solutions to the BLMP equation such as interacting solutions. For example, one can 
construct a rational-soliton (u rs ) and rational-positon (u rp ) solutions using Wronskian 
formulation [3J El El [11]. One gets 



and 



Ur,(-Z, t, 7l) 



lirp(z,*,7l) 



7i2tanh(7i(2; — 47ft)) — 1 

-27?* 
7i2;tan(7 1 (2; + 47 1 2 t)) + 1 



(20) 



(21) 





Figure 2. From left to right, a complexiton (|19[) (a = /? = 1), a rational-soliton (|20|) 
and a rational-positon (|2"Tj) solutions for z — x + q(y), t = 2 and 71 = 1. 

In the general case, one could be interested in solving the bilinear equation ([6j) and 
retrieve, for example, a general travelling wave solution. Indeed, if the following 
constraints are supposed 

1 



/ = f(w) = f (ax + ~ m (y) ~ ^aH^j > Uw = afv> + - « 2 )/- 



where a is an arbitrary complex constant, then the bilinear form ((6]) is automatically 
satisfied and we get a solution of the BLMP equation 



u 



-2ad w \og (c ie ^ 1+a)w + c 2 e 



m(y), 



(22) 
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where C\ and C2 are integration constants and m = m(y) is an arbitrary function of the 
variable y used in (jSJ) to bilinearize the BLMP equation . At our knowledge, the above 
solution is new. Here we give an example for a = and c\ — C2 = 1, one gets the 
kink-type solution 

w 

u(x, y, t) = — atanh(— ) — m(y). 

Most of the above solutions are new and will be useful in subsequent sections to 
construct superpartners of the supersymmetric BLMP equation. 



3. SUSY extension of the Bell polynomials 

The Bell polynomials are known as an alternative to tedious calculations of finding a 
bilinear form of an equation [Til 122] . Indeed, we are able to rewrite the equation in 
terms of the Bell polynomials and to get simultaneously the change of variable and the 
bilinear form of the studied equation. For more details, we refer the reader to [22J. 

For our present analysis, we introduce a particular case of the SUSY Bell 
polynomials. We refer the interested reader to [22] for a more general introduction 
to the SUSY Bell polynomials. 

Let A(X] 6) be a super different iable bosonic superfield [25] where X = (x,y,t). 
The super Bell polynomials are defined as 

n x , {kx , ky) (A) = Y, x , ikx , ky) (A, x , ikx , ky) ) = e~ A D^D k y yd^d l y y^e A ) (23) 

where I ■ X — (l x x, l y y, l t t) and D x and D y are the supercovariant derivatives defined in 

(ED. 

The super Bell polynomials enjoys a remarkable property given by 
Y (AM ._ Dt*D^d^f{X-9) 

yi-x,(k x ,k v ){A)\A=\nf - — - f — •- T(x~6) ' 

The super binary polynomials are defined, using the super Bell polynomials, as 

yi-X,(k x ,k y )(v,w) = Yi.x,(k x ,k y )(A s . X ,(n x ,n y )), (25) 

where the different derivatives of A are replaced by the superfields v and w using the 
following procedure 

' v s . x ,(n x ,n y ), s x + s y + s t + n x + n y is odd 
. w s . x ,(n x ,n y ), s x + s y + s t + n x + n y is even 



A s -X,(n x ,n v ) 



As an example, we have that 

^(3x,o,o),(o,o)(^) = A xxx + 3A X A XX + Al, 

3 ; (3a;,0,0),(0,0)(w, w) = V xxx + 3v x W xx + V^. 

The super binary Bell polynomials are related to the SUSY Hirota derivative 
Indeed, we have the following relation 

3U(^) (v = lnt w = l nf g\ = (fg)-iStS*yV l *V l y *V l *(f ■ g), (26) 
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where the super Hirota derivative is defined as 

S,(f-g) = (D,-D,,)f(X,9)g(X',9')\ x=x ,, e >= e . (27) 

The super P-polynomials |22] are a particular class of (121)]) where we take f = g. 
We have that 

Pi-x,(k x ,k y )M = yi.x,(k x ,k v )(0,w = 21ns) = g- 2 S k x *S k y yV l *V y yV\\g-g).{2%) 

A direct consequence of (|2T)j) is that Pixxk x ,k y )(w) = whenever / x + / y + l t + + k y is 
odd. 

4. Bilinear form, multisolitons and superpartners 

In this section, we propose a bilinear form of the SUSY BLMP equation ([7]) and give 
multisolitons and superpartner solutions. In order to get this form, let us cast the 
change of variable $ = D X F + A(y; 8) where F is a bosonic superfield, A is a fermionic 
superfield and integrate with respect to D x and secondly by D y . We thus get the 
fermionic equation 

D y (F t + DlF) - 3(D x F)(D y DlF) - 3(D y D x A)F xx = 0, (29) 

where the constants of integration as been set to zero and we recall that D x = d x . Using 
the following P-polynomials 

P(o,o,t),(o,i)(w) = D y w u -P(3x,o,o),(o,i)(w) = D y w xxx + 3w xx D y w x , 

equation ( 1291) can be rewritten as the polynomial equation 

-P(o,o,i),(o,i)(~-^) + -P(3:e,o,o),(o,i)(~-^) — 3(D y D x A)P^x,o,o),(o,o) (— F) = 0. 

We thus cast F = — 21ng and we get the bilinear form of the SUSY BLMP equation 
(J7J) given by 

(S y V t + SyVl - 3{D y D x A)Vl){g -g)=0. (30) 

Note that if A(y; 9) = rj(y) + 0c(y) then we get D y D x A = 6c' (y) and, in terms of 
components, if g = g b + Ogf, the bilinear equation fl3U]) becomes 

Qi{V x ,V y ,T>t){g h ■ g b ) = (V y (V t + Vl) - 3c'(y)V 2 x )(g b ■ g b ) = 0, 
Q 2 (V x ,V y ,V t )(g b ■ gf) = (V t + V 3 x )(g b ■ gt) = 0. 

In order to find the dispersion relations of g b and g* , one as to solve the set of equations 

Qi(k, p, u) = 0, Q 2 («, P, = 0, 
which is equivalent to 

p(u + k 3 ) - 3c'(y)n 2 = 0, u + k 3 = 0. 

Thus we get c'(y) = 0. This shows that to get non-trivial multisoliton solutions, one as 
to consider the super bilinear form 

S y (V t + Vl)(g-g) = 0. (31) 
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4-1. Multisolitons 

As usual, in order to get super multisolitons of the SUSY BLMP equation ((7|), we use 
the e-expansion of g 

oo 

g = 1 + 60i + e 2 g 2 + • • ■ = ^ e fc c/ fc 

fc=0 

and substitute this expansion in the super bilinear form (l3Tj) . For the super travelling 
wave or one super soliton solution, we take 

g 1 = e 01 , 0i = k x x + pxy + uit + 6& =<fi + 6& 

and Qi = for all i > 2. Introducing g in the super bilinear form ( 13TI) . we get the 
dispersive relation cji = — k\ . For the two super soliton solution, we take 

gx = e^ 1 + e 02 , g 2 = A 12 1 1 + 2 ^ 2 ) e iPi-Hpa+e(Ciai2-K2cwi) ( 32 ) 

V P2-P1J 

and gi = for all i > 3. In this case, we get the even relations for i — 1,2 

3 , _ (K1-K2XP1-P2) _ P» + Pj /ooN 
CJj — — Kj, A 12 — 7 ; 77 ; r, Olij — . [66) 

(«1 + ^2j(Pl + P2> Pi~ Pj 

The generalization to super multisoliton solutions is direct following the procedure 
described in [16] . An interesting feature of the above solution is the absence of fermionic 
constraints. 

We refer the readers to [121 [IS H3 US H3 [IS HE EH E31 IM] for more details on 
the construction of supersymmetric multisoliton solutions. 

4-2. Superpartners 

In section 2, we have presented numerous classical solutions of the BLMP equation (j2J). 
Here we get for some of them, the superpartner £ = (k, where ( is an odd constant and 
k = k(x, y, t) using $ = £ + 6u. 

Let us consider solution (122]) which takes the form 



/ \ 2cj(g?i — d 2 ) / \ / \ 

u [dl>d2) = -a(l + a) + di _ d2 + {di + d2)eW ~ m(yl (34) 

where d\ = c\ + c 2 and d 2 = C\ — c 2 . In this case, defining if) = £ w , we see that i/j = if)(w) 
satisfies the time independent Schrodinger equation [H, EJ [9] 

* (i m ~ dl) )i- (35) 

Vww ~ y 2(d a cosh(f ) + d 2 sinh(f )f ) [M) 

It can be solved using the change of variable \i = tanh(~). We get the solution £ = (k 
where ( 2 = and 

IV 3^/7 

k(di,d 2 )( w ) — /5i(4did 2 (lldi + dj) cosh — + lhdid 2 {d\ - d\) cosh — 
+5d 1 d 2 (dj - dj) cosh — + Adl(lldj + dl) sinh — 



-5(d? - d»)(4d? - 4) sinh ^ - (d? - d«)(4d? + 4) sinh ^)/ 

IV IV IV 

(di cosh — + d 2 sinh — ) + f3 2 (Adx{hd\ — d 2 )(<ii + ^2) cosn + 
5di(df — dl)(5d\ — 3d 2 ) cosh — + bd\(d\ — d 2 ) cosh — + 

IV 3?i7 

Ad 2 (5dj - d\){d\ + d 2 ) sinh — + 5d 2 (dj - t^) sinh — + 



d 2 (dj - d 2 2 )(9dl + d 2 ) sinh — )/{d x cosh — + d 2 sinh — ) + f3 3 , 



u {i,o) — — atanh(^) — aa — m(y) 



—)/{d 1 cosh - + d 2 smh — 

where the /Vs are constants of integration. 

Let us give some relevant special solutions. The case (di, d 2 ) = (1, 0) (corresponding 
to the Poschl- Teller potential |26j in the Schrddinger equation), we get the kink-type 
(soliton) solution 

and the corresponding superpartner = C^(i,o) where 

w 

k (1)0) (w) =fo- 8/3i(4sinh(w) + sinh(2w) - 2tanh(— )) 

+ 10/3 2 (4cosh(w) + cosh(2w)). 
The case (d\, d 2 ) = (0, 1) leads to 

u (o,i) — — acoth(— ) — aa — m{y) 
with corresponding superpartner £( ,i) = Cfyo.l) where 

fc(0,i)O) = (3 3 + 2(/3i - /3 2 )(-4cosh(w) + cosh(2w)). 

The above solutions have never been considered and we want to insist on the novelty 
which is due to a new approach of solving SUSY models. 

5. Bilinear Backhand transformations 

Backlund transformations PHElEimEOlEHEEJEghas been shown to be useful in 
constructing new solutions from known ones and give nonlinear superposition formulas. 

The above sections have exhibited new solutions of the classical and SUSY BLMP 
equation. Here we propose a nice application of the SUSY Bell polynomials applied 
to the SUSY case. Indeed, we propose to find bilinear Backlund transformations. To 
achieve this, we define 

^- / ( < ^ > ) D X Qyt ~\- D X ^ XXX y 3D X Dy(& xX DyQ X ~\~ ^ XX DyD X N^ 

and we suppose that ip = — 21nr and rj = — 2hiyU are such that L(ip) = L(rj) = 0. In 
order to use the super binary Bell polynomials, we define new variables 

„ = -*Z2 = h,(lV «, = -*+2 = l»( T/ ,) (36) 
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and, in terms of v and w, we have the following equation 

L(ip) - L{rj) = -2{D x v yt + D x v xxxy 
+3D x D y (v xx D y w x + w xx D y v x - v xx D y D x K)) = 0. 

In order to get bilinear Backlund transformations of the SUSY BLMP equation, we fix 
the suitable constraint 

y(x,o,o),(o,i)(v, w) - D y D x A = D y w x + v x D y v - D y D x A = 0. 
Then one can show, using the above constraint, that 

L(ip) - L{rj) = -2D x (y {0 ,o,t),(o,o)(v,w) + 3>(3*,o,o),(o,o) 
and thus we impose the polynomial constraints 

3^,0,0), (o,i) {v, w) - D y D x K = 0, y(o,o,*),(o,o) ( v , w ) + ^( 3a; ,o,o),(o,o) (v, w) = 0, 

which leads, together with the change of variables (1361) . to a set of bilinear transformation 

(S y V x - D y D x A)(r ■ n) = 0, (V t + V 3 x )(r ■ fi) = 0. (37) 

The above transformation is similar to the bilinear form of the SUSY mKdV equation 
proposed in [29] identifying the variable y with the variable x. It can be viewed 
as a SUSY generalization of a bilinear Miura transformation. Furthermore, taking 
r(x,y,t;8) = f(x,y,t) and fi(x,y,t;9) = g(x,y,t) in (!37j) . one gets the set of bilinear 
transformations obtained in [TT] for the classical BLMP equation given by 

(V y V x - c'(y))(f -g) = 0, (V t + V 3 x )(f ■ g) = 0. (38) 

We thus use this transformation to get a set of bilinear Backlund transformations 
of the SUSY BLMP equation. Indeed, we have the following proposition: 
Proposition Suppose that (r, fi) is a solution of the bilinear equations ( 1371) . then (r', //) 



satisfying the following bilinear relations 

V x {t ■ /i') - aV x {fi ■ t') = 0r/i' - apr'/i, (39) 

S y {r • ji) + aS y (fi ■ t') = jTfi' + (yyr'fi, (40) 

(V t + Vl- 3/3V 2 x + Zf3 2 V x ) (r • t') = 0, (41) 

(V t + V 3 X - 3/3V 2 x + 3(3 2 V X ) (ji ■ fi') = 0, (42) 



solves the bilinear equations (|37j) . where a and (3 are even constants and 7 is an odd 
constant. 

Proof: The proof of this proposition is similar to the one presented in [29]. Indeed, we 
take 

V x = 2 (T^'(S y V x - D y D x k)(r ■ fi) - Tfi(S y V x - D y D x A)(r f ■ fi')) , 
V 2 = r'fi'(V t + Vl)(r ■ fi) + Tfi(V t + V 3 x )(r' ■ fi'), 

and show that V\ = V2 = using the above bilinear relations. The proof follows noticing 
that V\ is equivalent to 

V x = 2 (r'fi'SyV x (T ■ fi) - TfiS y V x (r' •//))■ 
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The above Backlund transformations associated to the SUSY BLMP equation are 
related to the ones for the SUSY mKdV equation [29J and the M = 2 SUSY KdV 
equation [20] when one identifies the variable y with the variable x. 

6. Conclusion 

We have started this paper, with an extensive presentation of new solutions of the 
BLMP equation. Indeed, we have produced using Hirota and Wronskian methods 
numerous solutions, e.g rational, soliton, positon, negaton solutions. Most solutions 
being introduced for the first time. 

We also presented a new supersymmetric Boiti-Leon-Manna-Pempenelli equation. 
We have studied the SUSY BLMP equation from different view points using the beauty 
of the SUSY Bell polynomials. The novelty of this paper is the new approach consisting 
of finding the superpartners based on solving the linear equation ( 1331) 

Futur work should turn around the application of the super Bell polynomials to 
M = 2 extensions of nonlinear PDE's e.g. the SUSY KdV equation [30J. 
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